Intramolecular vibrational energy redistribution from a high frequency mode in the 
presence of an internal rotor: Classical thick-layer diffusion and quantum localization 
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We study the effect of an internal rotor on the classical and quantum intramolecular vibrational 
energy redistribution (IVR) dynamics of a model system with three degrees of freedom. The system 
is based on a Hamiltonian proposed by Martens and Reinhardt (J. Chem. Phys. 93, 5621 (1990).) 
to study IVR in the excited electronic state of para-fluorotoluene. We explicitly construct the state 
space and show, confirming the mechanism proposed by Martens and Reinhardt, that an excited high 
frequency mode relaxes via diffusion along a thick layer of chaos created by the low frequency-rotor 
interactions. However, the corresponding quantum dynamics exhibits no appreciable relaxation of 
the high frequency mode. We attribute the quantum suppression of the classical thick-layer diffusion 
to the rotor selection rules and, possibly, dynamical localization effects. 
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Intramolecular vibrational energy redistribution (IVR) 
has been the focus of several theoretical and experi- 
mental studies for the past few decadesii^i^i^i^. An 
important challenge, from the viewpoint of prediction 
and control of reaction dynamics, is to understand 
the influence of specific structural features (functional 
groups) of a molecule on the IVR dynamics. Several 
studies^'^'^'^ i-^^'^^i^^'^'^i-^'^i-^^ suggest that the early stages 
of IVR can indeed be described from such a moiety- 
specific i.e., local viewpoint. In particular, starting with 
the experiments by Parmenter and coworkersiSiiiii^iii^ 
the role of large amplitude motions involving internal ro- 
tors has received^>i^'i^>^i^>^'2^i^'^>2&'^, and continues 
to receive^i^'^i^'^, considerable attention in the liter- 
ature. 

In this work we focus on a mechanism, suggested by 
Martens and Reinhardi*2& and inspired by the observa- 
tions of Parmenter and coworker a^^'^^ , for the accelarated 
IVR rates among the ring modes in para-fluorotoluene 
(pFT) in comparison to a similar molecule without the 
methyl rotor (para-difluorobenzene). The essential fea- 
tures of the proposed mechanism^^, based on detailed 
classical dynamical studies of a model Hamiltonian for 
pFT incorporating the methyl rotor, five of the lowest 
frequency ring normal modes, and important ring- rotor 
couplings, can be illustrated using the 'minimal' Hamil- 
tonian 
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In the above, three degrees of freedom Hamiltonian, the 
ring normal modes 1, and 2 (wi < uj^) are coupled with 
the hindered rotor mode but the oscillators are not di- 
rectly coupled to each other. Setting ai(a2) = one 
obtains two degrees of freedom subsystem describing the 
high (low) frequency mode-rotor coupling. The subsys- 
tem dynamics can be analyzed as follows. A resonance 
LOj = iujr, driven by the coupling terms in Eq. ([T]), occurs 
at Pp' « ±LUj/6B and leads to quasiperiodic exchange of 



energy between the oscillator and the rotor. However, if 
in addition the condition pjjj'^ ~ -y/Vb/i? is satisfied then 
the LUj = 3ujr resonance can overlap with the hindered 
rotor separatrix, Er — Vq, generating large scale stochas- 
ticity in the phase space leading to a rapid relaxation of 
the oscillator energy. Thus one anticipates that if at a 
given energy the low frequency-rotor subsystem (q!2 — 0, 
"bath") has a strongly chaotic phase space then an ex- 
cited low frequency mode will relax rapidly. On the other 
hand at the same energy the high frequency-rotor sub- 
system (ai = 0), typically, would have a near-integrable 
phase space and hence an excited high frequency mode 
will not relax. 
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FIG. 1: Illustrating the mechanism for relaxation of a high 
frequency mode (2) according to Martens and Reinhardt^. 
In (A) the frequency loi is sufficiently high such that the m:l 
resonance between the bath (1) and the rotor (r) does not 
overlap with the hindered rotor separatrix Er = Vq. On the 
other hand, in case (B) the bath frequency is low enough for 
an overlap resulting in large scale stochasticity. This band 
of chaos now leads to the relaxation of the high frequency 
system over long time scales. 

Interestingly, Martens and Reinhardt observed^^ that 
coupling the subsystems i.e., ai,a2 ^ leads to a slow 
relaxation of the excited high frequency mode despite the 
lack of resonance overlap. It was suggested^^ that this 
could be due to a slow diffusion along the thick layer 
of chaos created by the bath subsystem (cf. Fig. [TJ3). 
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Evidence for such a "stochastic pumping" mechanism^^ 
was provided by modifying the system so that the bath 
dynamics is no longer strongly chaotic as sketched in 
Fig. [T]A_. Clearly, such a mechanism would lead to exten- 
sive mixing of all the ring modes and therefore support 
the experimental observations. However, there are cru- 
cial questions, given the debate about the precise role of 
the rotor— that have not yet been answered. Firstly, 
the Hamiltonian used in the original studies involved five 
of the ring modes and the hindered rotor. Is it possible 
that the 'minimal' Hamiltonian in Eq. ([T]) is enough to 
capture the main mechanisms? Secondly, can one explic- 
itly show the existence of the thick-layer diffusion and 
understand the nature of this diffusion? Finally, is the 
mechanism valid quantum mechanically? The last issue 
is crucial, as was also hinted by Martens and Reinhardt, 
due to strict rotor selection rule o^^i^^ and the possibil- 
ity that quantum mechanics might localize the long time 
classical transport. In this work we report some of our 
preliminary results which provide answers to the ques- 
tions posed above. 

We consider the Hamiltonian in Eq. ([T]) with the os- 
cillator frequencies lui = 110 cm"-'^, 0^2 = 359 cm~^, ro- 
tor barrier Vq = 34 cm~^, and methyl rotational con- 
stant B ~ 4.65 cm^^. The parameters are identical to 
that of the earlier worl^ except that only the lowest 
and the highest frequency modes have been retained'^^. 
An ensemble of classical trajectories with fixed actions 
{Ji, J2, P4,) and random phases of the oscillators (6'i, 6*2) 
and the rotor (0) are generated to model initially excited 
nonstationary zeroth-order quantum states \vi,V2,Vr) ■ 
In Fig. [2] we show the results for a representative case 
of the excited high frequency oscillator corresponding to 
the zeroth-order quantum state |0, 5, 1). Trajectories are 
propagated for 100 ps and in Fig. [5] the ensemble aver- 
aged mode energies are shown as a function of time. It 
is clear that the high frequency mode relaxes over the 
long timescale with the energy being gained by the low 
frequency mode and the rotor. In fact the subsystem 
rotor surface of sections shown in Fig. [2] correspond to 
H = i?05i ~ 2034 cm~^ and suggest that the high fre- 
quency mode should not relax. Thus the energy flow out 
of mode 2 in Fig. [2^ is a three degree of freedom effect 
and identical to the observations made by Martens and 
Reinhardt in their six degrees of freedom Hamiltonian^^. 

Is the long time energy relaxation of the mode 2 seen 
in Fig. [2^ due to the thick-layer diffusion mechanism sug- 
gested earlier^^? To this end we study the classical IVR 
dynamics by explicitly constructing the zeroth-order ac- 
tion space {Ji, J2, Pip) state space. A set of initial 
actions were chosen by fixing the total energy H — Eq^i 
and angles {9i, 62, 4>) = (0, 0, 0). These set of initial con- 
ditions are then propagated to T = 100 ps and the diffu- 
sion constants 

dkiT) = ^ [ \nkit)~n\dt (2) 

^ Jo 

for = 1, 2, r are determined. The Qk are the local non- 
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FIG. 2: (Color online) (a) Classical average mode energies 
(in cm^^) versus time for an ensemble of trajectories with 
initial conditions corresponding to \vi,V2,Vr) — |0, 5, 1). The 
coupling constants^^ are ai = —3.14 cm~^/A and 02 ~ 2.3 
cm~^/A. The (0, P^) surface of sections, at the zeroth-order 
energy -B051 = 2034 cm~^, for the two degrees of freedom 
subsystems obtained by decoupling the high and the low fre- 
quency modes respectively are also shown. The 3:1 resonance 
islands (green) and the rotor (blue) regions are clearly visible 
in the 2-r subsystem, (b) The (6I1, 6*2,0) = (0,0,0) slice of 
the state space at £051 is shown and exhibits regions with 
large and small diffusions (in ps~^). Note the existence of 
a chaotic band which should be compared to the schematic 
shown in Fig. [1)3. A representative trajectory with large dif- 
fusion along the chaotic band is also shown in (b) with the 
circles representing 1 ps intervals, (c) The frequency content 
of the representative trajectory in (b) confirms the chaotic 
behaviour of the rotor-low frequency subsystem. See text for 
details. 



linear frequencies extracted using a wavelet-based time- 
frequency approach^. The sum of the d^s is then asso- 
ciated with each trajectory (point in the state space) and 
hence the regular (low diffusion) and chaotic (high dif- 
fusion) regions are identified. In Fig. [^b the state space 
is shown as a projection onto the (Ji, J2) plane and the 
band of stochasticity, due to the mode 1 and rotor in- 
teraction as sketched in Fig. [1)3, is clearly seen. Several 
classical trajectories corresponding to |0, 5, 1) were prop- 
agated for 100 ps and a variety of diffusive behaviours 
were seen. In Fig. ^jp we show a typical high diffu- 
sion trajectory superimposed on the state space which 
clearly shows the diffusion along the thick layer of chaos. 
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This establishes the mechanism proposed by Martens and 
Reinhardt for the minimal Hamiltonian in Eq. (jT]) . More- 
over, to the best of our knowledge, this is the first time 
that the thick-layer diffusion has been explicitly shown 
in a molecular system. 
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FIG. 3: (Color online) Quantum average mode energes (red) 
compared to the classical (black symbols) analogs shown in 
Fig. [2^. The bottom panels compare the classical (thick 
black) and quantum (red) rotor angular momentum distri- 
butions, (a) Original system (cf. Eq.[l| indicates localization 
of the quantum p{P,p) as compared to the classical case. In 
(b) and (c) weak couplings between mode 1 and the rotor of 
the form f3Qism2(j> and PQisincf) with (3 — —1 cm~^/A re- 
spectively are added to the original system in order to break 
the Avr = ±3 selection rule in (a). 

Certain features of the dynamics that are inherently 
due to the system having three or more degrees of free- 
dom are worth pointing out at this stage. A careful look 
at the trajectory shown in Fig. ^jp reveals that for nearly 
15 ps there is very little energy relaxation in £2- In our 
computations this feature was common to a large num- 
ber of trajectories. A preliminary understanding of the 

15 ps 'incubation' time can be gained by the time- 
frequency analysis^^ of the classical dynamics. In Fig. [2]: 
the results for the representative trajectory undergoing 
thick-layer diffusion in the state space are shown as spec- 
trograms i.e., the distribution of the local frequency am- 
plitudes as a function of time. One can see that during 
the incubation time the low frequency mode and the ro- 
tor exhibit strong coupling and chaotic dynamics while 
the high frequency mode is decoupled. At longer times 



perturbations in the local frequency associated with the 
mode 2 are clearly visible. The rotor phase space shown 
in Fig.[2]for the bath indicates a thin layer of chaos which 
is insuffcient to drive the energy relaxation of the high 
frequency mode. Consequently we interpret the incuba- 
tion time as the time required for the low frequency-rotor 
subsystem to generate the thick layer of chaos. Inciden- 
tally, we remark here that Martens and Reinhardt arrive 
at nearly the same timescale for the high-frequency cut- 
off in the spectral density associated with the random 
force generated by the chaotic bath dynamics. After the 
incubation time the trajectory rapidly moves along the 
chaotic layer resulting in the energy flowing into the low 
frequency mode and the rotor. At longer times the trajec- 
tory once again performs diffusive motion about specific 
regions in the state space. 

We now address the crucial issue of wether the quan- 
tum dynamics of Eq. ^ would also exhibit the thick- 
layer diffusion leading to the relaxation of the excited 
high frequency mode. Surprisingly, and as far as we can 
tell, there has been no attempt to confirm the mecha- 
nism proposed by Martens and Reinhardt in the quantum 
domain. We note, however, that Leitner and Wolynes 
have studied22,!^ the effect of quantization on the weak 
stochastic pump model of Arnol'd diffusion. In this case, 
involving anharmonic oscillators, it was found that the 
quantization limits the extent of diffusion. Does this im- 
ply a similar localization for the quantum dynamics of 
our minimal Hamiltonian? To investigate this possibility 
we have performed detailed quantum studies of the IVR 
dynamics of our system. A converged basis set, using 
direct product of harmonic oscillator and free rotor basis 
functions, was employed in all the studies. In Fig. [3] we 
show the quantum average mode energies (Ek) for the 
zeroth-order state |0, 5, 1) and compare with the classical 
results of Fig. [51 It is clear from the figure that quan- 
tum mechanically the high frequency mode does not relax 
even on the 100 ps timescale - infact a beating pattern is 
observed. It is important to mention that results of fur- 
ther studies (not shown here) on the quantum dynamics 
including the rest of the ring modes, as in the original 
study^^ , again indicate little to no energy flow out of the 
high frequency mode. Thus, it appears that the classi- 
cal thick-layer diffusion is completely suppressed in the 
corresponding quantum system. 

There can be a couple of reasons for the observed quan- 
tum behaviour. The first, as already alluded to earlier, 
could be due to a genuine quantum localization of the 
thick-layer diffusion. Consequently, relevant eigenstates 
of Eq. ([T]) were investigated and none of them displayed 
a strong delocalization along the chaotic layer. Never- 
theless, the fact that the high frequency mode does not 
relax even upon the inclusion of additional ring modes 
points to an insufficient density of states in the chaotic 
layer. One other possibility has to do with the fact that 
quantum mechanically there is a strict rotor selection rule 
Avr = ±3 due to the Q sin 30 coupling in Eq. ([T]) . Thus 
the low frequency-rotor coupling cannot lead to the pop- 
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ulation of specific rotor states as opposed to the situa- 
tion in the classical dynamics. The consequences can be 
clearly seen in Fig. [3K where the quantum rotor momen- 
tum distribution^^ 

p{Vr = hP^)^^ I dtY,\{viV2Vr\e-'"''''\Q,b,l)'^ 

(3) 

is compared with the corresponding classical 
distribution^^ p{P^) at T = 100 ps. The quantum 
distribution is localized and the asymmetric form of 
the envelope is due to the selection rule. In order to 
assess the importance of the rotor selection rule to the 
observed quantum localization weak coupling terms 
between the low frequency mode and the rotor of the 
form l3Qi sin 20 and PQi sin were introduced in Eq. ^ 
with (3 = —1 cm~^/A. The resulting quantum and 
classical ^(P^) are compared in Fig. and Fig. [3J; 
respectively. Although one observes population of higher 



momentum states, the quantum distribution is still 
narrower than the classical distribution. At the same 
time even for the modified cases no energy relaxation 
of the high frequency oscillator is observed. Therefore, 
given the lack of a sufficient density of states in the 
chaotic layer, rotor selection rule might not be the sole 
factor responsible for the quantum localization. The 
classical-quantum correspondence of p{P(t,) have been 
studied earlier— in the context of periodically kicked 
one degree of freedom systems. In our case the system 
has three degrees of freedom and the nature of partial 
barriers is not yet clear. Further work, in addition to a 
careful study of the dependence of energy diffusion and 
quantum eigenstates on ?i, is required to provide a firm 
answer and such studies are currently underway. 

The authors acknowledge Profs. Gerrit Groenenboom, 
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